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ABSTRACT 

Motivated by the possibility of radiation driven instabilities in rotating magnetic stars, 
we study the stability properties of general linear perturbations of a stationary and 
axisymmetric, infinitely conducting perfect fluid configuration threaded by a magnetic 
field and surrounded by vacuum. We develop a Lagrangian perturbation framework 
which enables us to formulate a strict stability criterion based on the notion of a 
canonical energy (a functional of the fluid displacement ^ and its first time deriva- 
tive). For any given choice of {C, 9t^}, the sign of the canonical energy determines 
whether the configuration is stable or not at the linear level. Our analysis provides 
the first complete description of the stability problem for a magnetic star, allowing 
for both rotation and the presence of a magnetic field in the exterior vacuum region. 
A key feature of the Lagrangian formulation is the existence of so-called 'trivial' fiuid 
displacements, which do not represent true physical perturbations. In order for the 
stability criterion to make rigorous sense one has to isolate these trivials and consider 
only the physical 'canonical' displacements. We discuss this problem and formulate a 
condition which must be satisfied by all canonical displacements. Having obtained a 
well-defined stability criterion we provide examples which indicate that the magnetic 
field has a stabilising effect on radiation driven instabilities. 



1 INTRODUCTION 

The stability of a hydromagnetic configuration under the influence of linear perturbations is a classic interdisciplinary problem 
who se relev a nce ra nges from astrophysical systems to terrestrial laboratories. Accordingly, the literature on the subject is vast 
(see iMestell ^20oi ) for an up-to-date discussion of magnetic astrophysical systems and a rich list of relevant references). The 
present paper is motivated by the interesting possibility that the emission of radiation may drive various stellar oscillation 
modes unstable. This class of instabilitie s has so far only been studied in detail for systems radiating gravitational waves at 
a significant level (see lAnderssonI \20oi ) for a review). A very important question concerns to what extent electromagnetic 
radiation is relevant in this context. To begin addressing this question, we focus on the stability of astrophysical hydromagnetic 
systems, self-gravitating fiuid bodies endowed with both rotation and magnetic fields. Neutron stars are prime examples of 
such systems, as they are the most rapidly spinning and highly magnetised systems known in the cosmos. 

The found ations for the modern stability theory for Newtonian fluid stars were laid by Chandrasekh ar more than 50 



years ago (see IChandrasekhaJ (|l96ll ) for a summary). A major contribution to the subject was made by iBernstein et al 



( 1958ll who devised a stability criterion, based on the sign of the potential energy, for a static hydro magnetic configuration 
where the unperturbed magnetic field is tangential at the fluid- vacuum interface. Subsequent work bv lFrieman fc Rotenberj 
generalised the analysis to stationary conflgurations (albeit only for flelds that are conflned to the fluid). Another key 
generalisation was provided by Kovet j ( 19661 ) who considered the static problem but allowed the magnetic field to extend into 
the vacuum domain without any restrictions. 



For fluid configurations, an important contribution was made bv lLvnden-Bell &: Ostrikeij l|l967f ). Their stability criterion 
for a rotating fluid conflguration was expressed in terms of an operator which depends on the L agrangian displacemen t 
^. Since this oper ator represents the system's potential energy, this criterion is identical to that of Bernstein et al. ( 19581 ). 
A few years later Friedman fc Schutz ( 1978ah added a crucial element that had been overlooked in the previous analyses. 
They reformulated the stability criterion in terms of the Lagrangian canonical energy functional Ec{^,dt^), emphasising the 
existence of 'trivial' fluid displacements. Even though these triv ial displacements produce n o physical change in the system, 
they do change the canonical energy. Consequently, as argued by Friedman fc Schut j (1978a), the presence of trivials disables 
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the stability criterion of Lynden-Bell & Ostriker. To avoid this impasse iFriedman fc Schuta l|l978af ) identified the family of 
trivial displacements, and expressed the stability criterion in terms of physical 'canonical' displacements. Since the late 1970s 
this st ability criterion for rotating fluid configurations has been the final word on the topic. A particularly important outcome 
of the Friedman fc Schut j 1 1978al ) work was the discovery of the so-called CPS instability by which a l l rotating perfec t fluid 
stars are secularly unstable due to the emission of gravitational radiation ( Friedman fc Schutz 1978b : Andersson! 20031 ) . 

Quite surprisingly, to the b est of our knowledge , there has so far bee n no attempt to modi fy the stability criteria for 
magnetic sta rs along the l ines of Friedman fc Schutz ( 1978a[ l. The original Bernstein et al. ( 1958h criterion is now textbook 
material (see Mestel 1 2003h ). Yet it is clear that the existence of trivial displacements may invalidate the analysis for rotating 
stars. In fact, since they are characterised by a non-vanishing canonical energy, trivial perturbations (which are totally 
unphysical) might well be qualified as being unstable. This is troublesome, and highlights the need for the problem to be 
revisited. 

In this paper we aim to provide a rigorous stability criterion for rotating hydromagnetic configurations taking due 
account of the presence of trivial displacements. We consider a self-gravitating perfect fiuid threaded by a magnetic field, 
assuming a stationary equilibrium and an arbitrary st ructure for the magnetic field (extendable to the vacuum exte rior). 
This set up combines and generalises previous studies ( Bernstein et al. 19581 : Frieman fc Rotenberg 1960l : Kovet j 19661 ) . We 
determine the conditions that govern the family of trivial displacements and provide explicit conditions that need to be satisfied 
by physically acceptable canonical perturbations. The required canonical energy functional is also explicitly formulated. 
Given these results, the hydromagnetic stability criterion can be properly formulated and evaluated for any magnetic field 
configuration. To illustrate this, we apply the stability criterion to a few simple cases. In particular, we demonstrate that 
there are always canonical data which make a magnetic rotating star secularly unstable, under the emission of electromagnetic 
and/or gravitational radiation, in the limit of m ^ 1 (assuming that the perturbations depend on the azimuthal angle as 



2 LAGRANGIAN PERTURBATION THEORY FOR MAGNETIC STARS 
2.1 Formulation 

We consider a stationary, but otherwise arbitrary, self-gravitating perfect fiuid configuration (of compact support, surrounded 
by vacuum) which is threaded by a magnetic field B". The fluid flow is described by the velocity field u" and furthermore, ideal 
MHD conditions (i.e. inflnite conductivity) and a one-parameter equation of state p — p{p) are assumed. We will carry out 
the analysis in a coordinate basis, distinguishing between co- and contra-variant objects. This may be a somewhat u nfamiliar 
appro ach to a "standard" Newtonian fluid dynamics problem, but it has significant advantages (see lAndersson fc Comer 



1 20061 ) for a recent discussion). Most importantly, it allows us to define Lagrangian perturbations and discuss conservation 



laws in terms of the Lie derivative. The advantage of this will be clear at various points of our discussion, eg. when we consider 
the conditions at the surface of the star. As usual in this context spatial indices are raised and lowered with the fiat metric 
Qab, i.e. we have Ua ~ Qabu'' etcetera. 

The coupled fluid-magnetic field system is governed by the familiar Euler and induction equations (denoting the convective 
derivative by D/Dt = dt + u*Vb), 

^ + -Va(p+^]+S/a^~-^B'\/tBa^O (1) 

Dt p \ STV J i-Kp 

{B''Vt)Ua~Ba{Vbu) (2) 



Dt 

together with the continuity equation, 

dtp + Vaipu'') =0 (3) 
and Poisson's equation for the gravitational potential, 

= 4tvGp (4) 
Additional equations are Ohm's law (which determines the electric field within the fiuid), 

= -ie'""=M6Bc (5) 
c 

and, of course, the fact that VaS" = 0. 

The stationary 'background' configuration is obtained from the above system by setting all time derivatives to zero. We 
wish to consider an arbitrary perturbation of this background configuration. For the perturbative quantities we choose to 
work in the Lagrangian sense, by means of a displacement vector By definition we know that. 
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Au" = dtC 

For the covariant component we have, 

AUa = dt^a + u'Vt^a + u''V a^b 

where we have used 



For an arbitrary tensor Q (suppressing indices) we have the general expression, (jPriedman fc Schutall978al l. 
AQ = SQ + CfQ 

relating Lagrangian and Eulerian perturbations. For a scalar quantity p the Lie derivative is given by 
Meanwhile, we have 

^^Ua = S!'VbUa + UbVaS,'' 

and 

L^u — 5 VbU — u Vbc, 



(6) 
(7) 
(8) 
(9) 

(10) 

(11) 
(12) 



for CO- and contravariant vectors. 

Our first task is to express Ap , Ap, Ai?" in t erms of the displacement From conservation of mass and the equation 
of state we readily find (jFriedman fc Schutall978al l. 



Ap = -pv.r, 



Ap = ^Ap 



(13) 



where 7 — {p/p){dp/dp) is the adiabatic index. 

In order to find AB" we can perturb the induction equation ([2]). The procedure becomes more elegant if we first rewrite 
this equation (with the help of eqn. ([3])) using the Lie-derivative, 



[dt + Cn] 



P 







(14) 



This relation shows that B'^ / p is preserved along the fiow This encodes the expectation that the field lines are 'frozen 
into' the fluid. An attractive feature of (|14[) is that the A and [dt + Cu\ operators commute. Hence, we immediately have 



[dt + Cu]A 



P 







(15) 



In other words, the Lagrangian perturbation of B°' j p is also Lie-transported by the flow. It is natural to assume that the 
initial conflguration is unperturbed, i.e. to choose the trivial solutior0 A(_B"/p) = to (|15p . This then leads to 



AB' = -fi'Vcf 



From this we readily obtain, 

ABa = AgabB'' + QabAB'' 



ABa = BbVat - BaVb^" + B''\/ b^a 



The corresponding Eulerian perturbation of the magnetic field is given by 
b" = SB" = AB" ~ C^B" ^ b" ^ -B'^iVb^'') - (C'V^B" + {B''\/b)t 
This is equivalent to (jRobertsI 1967 ). 

where e"'"^ is the familiar three-dimensional Levi-Civita symbol. Note the useful identity, 

abc cb fC c-b r-c 

e eade = OdOe ^ <JeOd 



(16) 
(17) 
(18) 
(19) 
(20) 



^ It is worth noting that this choice ensures that the Eulerian perturb ation of the induction equat ion ||2j vanishes automatically, as 
required for an unconstrained variational formalism, see the discussion in [Friedman fc Schutj l ll978al'l . 
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2.2 Master equation 



We now have all the tools required to derive a 'master' perturbation equation involving only the displacement vector 
This can be achieved by considering the Lagrangian perturbation of the Euler equation ^ . With the help of the background 
equations and after some rearrangement we end up with 



pdUa + 2pu''Vbdt^a + p{uVbf^a + 5 p{ V + uVcUa ) + pVa5$ + V. (Sp + B, 

\ 47r 



4:71 

This equation is of the form, 



(21) 



(22) 



where A,B,C are operators constructed fr om background variables. T he first two are identical to the operators appearing 
in the non-magnetic limit and, as shown bv lFriedman fc Schuta (jl978al ). they are respectively Hermitian and anti-Hermitian 
with respect to the inner product, 

(,7^^a> = jinT^adv (23) 

The integral is taken over the space occupied by the fluid and the asterisk denotes complex a conjugate. That is, we have 
(r;,A[e]) = (C,^W)*, (^,B[e]> = -(e,5W>* (24) 



where both r;" and are solutions of eqn. (I22|l . Following [Friedman fc Schuta (|l978al ) we can construct the 'symplectic' 
product, 

W/(r;,C) = {77,n(0>-(n(77),O (25) 
where we have introduced the conjugate canonical momentum, 

m) = mi\ + \Bm = p^ (26) 

Then it follows (using eqns. ([22]) and ([24])) that, 

dtW{7^,0 = {CM,0-{ri,cm} (27) 

This shows that, provided that t he operator C is Hermitian, W{ri,^) is a conserved quantity. This is true for the non-magnetic 



case IjFriedman fc Schutall978al ). and we can show that it remains true for the magnetic case as well. 

The proof of the Hermiticity of the C operator is somewhat lengthy. In order to avoid overloading the reader with algebra 
at this point, we relegate this calculation to Appendix|X] The final result is, cf. Eq. (|A15|I . 



r7"(5p + B"bc/47r) 



47r 



V47rG 



-^V'^^^Vft^^* + ^B'B"Ver?"V6ea + v't'^aVtip + B^/8n) + V^ri^VbthP + B^/A-k] 
+ 4^ [V.77'^B'=(e'VcB5 - VcC'Sb) + V6C''B'('7''VcBa - Vc7?''Ba)] 

-p[uV,Ua + V,<E>)(r?"V6e' + r Vi,??') + py^^C'V, Vi,$ - pUu''ViT)''SJ,^a 



(28) 



The fact that rj"" and ^" appear in a symmetric way proves that the operator C is Hermitian, modulo the surface terms 
(resulting from the total divergence in the first line of (|28p '). To discuss these remaining terms we write the result as, 

(29) 



where Ch is Hermitian and the surface terms are given by. 



Sm(^7,0 



V 47rG 



(30) 



(31) 



In the _B — > limit eqn. (|29p reduces to the result of iFriedman fc Schuta (|l978al \ In this case the fluid surface terms vanish 
identically since we have p = p = Vp = and Ap = at the sur face. The remaining gravitational term g'^''S^$\/b5n$ can be 
shown to be Hermitian, see Andersson. Comer fc Grosart 1 2004h . rendering the full C Hermitian. Using the same arguments 
it follows that the Sf term is Hermitian also when the magnetic field is present. The remaining surface term is, however, 
generically non zero. This is easy to appreciate since the magnetic pressure does not have to vanish at the fluid surface. A 
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simple way to avoid problems would be to assume that the surface terms vanish anyway, e.g. by only considering magnetic 
fields that are confined to the fluid. This is, however, a severe and not very realistic restriction. Instead, we will consider the 
general problem and prove that the remaining surface terms are Hermitian. This requires further manipulation and a detailed 
discussion of the appropriate boundary conditions that need to be satisfied at the fluid/vacuum interface. 



2.3 Surface boundary conditions 

Since we are interested in studying general systems which have a fluid- vacuum interface, we need to consider the relevant 
boundary conditions at such a surface. In order to work out these conditions we first consider the normal vector to the surface 
and how it is affected by a Lagrangian perturbation. The Lie d erivative plays a key role in this analysis and simplifies the 
problem considerably, (compared to, for example, the analysis in iRobertd ([1967|)). Let us identify the surface as a level set of 
a scalar function /. In the fluid it would be natural to identify this function with the fluid pressure, but we also assume that 
the function can be extended to the exterior domain in a smooth fashion (this extension facilitates the deflnition of derivatives 
of /). We then require, 

\dt + £u] / = (32) 
from which it is easy to show that, 

[dt+ Cu]Vai = Q (33) 
In other words, the gradient Va/ is constant in the fluid frame. We also see that we must have, 

\dt + £u]A/ = (34) 

The trivial solution to this equation is A/ = 0, which essentially means that a fluid element at the original surface remains 
at the perturbed surface. The normal to the surface can obviously be taken to be = ^af, from which a unit normal vector 
can be constructed, 

n. = ^, where iV = (g^^V^/V,/)^/^ (35) 
Given the above results, it is straightforward to show that 

[dt + jCu]na = -^[dt + Cu\N (36) 

This shows that any change to the unit normal is parallel to the normal itself. Hence the orthogonal projection (into the 
surface) is conserved by the flow. This leads to, 

= (n n - g )n VbUc (37) 

which is the equation for the normal vector given by, for example, Kovet j ( 19661 ). Perturbing this equation leads to. 

Aria = hahcnVbi" (38) 

In deriving this we have assumed that the background configuration is such that h°'Ua ~ 0, that is, there is no local expansion 
or contraction. Equation (|38[) highlights a key difference between the normal vectors ha and Ua, since for the latter we have 
from 



[dt + Cu]Ana = (39) 
It is natural to use the trivial solution 

Ana = (40) 
for fluid elements on the surface. Then it immediately follows that 

g^'^'Aui = An" + n^V^C + ff" V = (41) 
That is. 

An" = -n^VcC - g'^nVi (42) 

Having a mathematical representation of the perturbed surface at hand, we can move on to the surfa ce condit i on for the 
magnetic and electric flelds. The general conditions can be found in standard textbooks, see for example. [Roberts! (119671 ). In 
describing them we will use angular brackets to represent the step in a quantity as we pass from the fluid to the exterior. 
That is, for any quantity Q we define, 

(Q) = Q^~Q (43) 
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where the index 'x' designates an exterior vacuum quantity. The brackets denoting a jump across the boundary should not be 
confused with the brackets used in the definition of the inner product (|23p . One can easily avoid mixing things up by keeping 
in mind that the symplectic product takes two arguments. 

From the fact that B" is divergence-free we immediately obtain, 

n''{Ba)=0 (44) 

Meanwhile, the Maxwell equation e'''""VbEc = ~{l/c)dtB'' leads to, 

e^'^ntiE,) = -(h''ut){B'') (45) 



The surface conditions (|44|l and (I45p are originally formulated in terms of the unit normal vector n". However in view of 
the previous results (eqns. (I38|l and (|40p ) the forthcoming analysis is greatly simplified if we rewrit^ the surface conditions 
in terms of = Nh'^. Subsequently taking the Lagrangian variation of (|45p . and using the fact that Uau"' = for the 
background, as well as eqn. (|4U|I . we find 

ch,A{e''''^E,))^hidtt{B'') (46) 
Use of the background equation ([5]) for the interior electric field together with eqn. (|16|l allows us to show that 

Aie"'"' Ec) = - i [{dtCB'' - dt^B") + [v!' B" - u''B'')\/cC] (47) 
For the exterior electric field we have 

A(e"''"£j) = -e^'^B^VfeC* + e^'^AB? (48) 



where we have used ( Friedman fc Schut j 1978a ) 



Given these relations, eqn. (|46p becomes 

ce-^'^n^AEl = nfc(ac'Bx - B'^dtC) + uUbB^VuS!' (50) 
where we have used (for the background) 

e"'"=n6-B: = (51) 
which follows from (I45II. 



Finally, it is useful to express the exterior electric field in terms of the vector potential , 
El = --^dtAk (52) 
Perturbing this we find 

AEl = --dt{AAk - C^Ak) = --dtSAk = --dtAk (53) 
c c c 

where we have used the fact that the background is stationary. This leads to the final form for the required boundary condition 
nbdtCB'l, ~ nbB'dtC + ^"""utdtAc = (54) 
This relation is a total time derivative which upon integration gives 

nt^Bi ~ h.B'C + e'""'h,Ac = (55) 
This is identical to the condition derived by Kovet j (196^) in the case of a static background configuration. 



An additional interface condition that needs to be satisfied is the continuity of the normal component of the fiuid-magnetic 
field stress tensor, 

( T,t,n ) = 0, where = -g^i ( P + ^ ) + 4^ ^6 (56) 



Explicitly, 
E 

Btt / 47r ^ 



"-(p+lr ) = i(^^5')(S'"> (57) 



Another way of seeing that the perturbed norm of the normal vector does not affect the result is to note that, when we perturb l|45|l 
then AA'' will multiply a quantity that vanishes since the boundary conditions also hold for the background. 
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Projecting this along n" and using (I44|) we obtain, 



P + 



8n 



= 



and 



{nBt){Ba)=0 

Lagrangian variation of these conditions (replacing with n" in (|59|) . as before) results in, 

Ap = ^ ( + BJC' VtB^ - B.fe^ + B^'VtB' ) 

and 

ncB'(b'"-bS-C''V5(B'")) =0 



(58) 



(59) 



(60) 



(61) 



It is worth emphasising that we are assuming that ideal MHD can be used throughout our system. This assumption will 
not be valid for a realistic neutron star model since the Alfven velocity will diverge near the fluid surface, as p — > 0. A detailed 
analysis of the surface region would require the introduction of a transition region from ideal MHD conditions in the core of 
the star to the full Maxwell equations in the exterior vacuum. Such an analysis is, however, much more intricate than our 
slightly inconsistent approach. Once the stability problem has been completely understood within the present framework, it 
would be interesting to extend the analysis to a more detailed surface/vacuum model. 



2.4 Magnetic surface term 



Returning to eqn. (|29|) . our immediate objective is to e liminate , or sy mmetrise Sm, thus rendering the operator C Hermitian. 
The following steps closely resemble the calculation bv iKovetz j 19661 ). 
For convenience we split Sm into two pieces: V ■ Sm ~ Si + S2 with. 



Si 



S2 = Va [B'"r?"(bc + C'V6Bc)] 

We write Sp — Ap — ^''VbP, which with the help of eqns. (I60|) and (|62|) leads to 

Sl=Va 



Since eqn. (|58p is uniformly valid everywhere at the fluid- vacuum interface we have 
B^ 



ahc T—7 I I 
e 726 Vc ( P + 







The cross product of this vector with is also zero, which leads to 
n.t^,i^P^^^^h,tsj^^p+^ 
Hence, 

h(ri,^) = J dVSi = J dS ^anbrj''^''h''Vc{p + ^ ) + -^fiar]'' B^b"^ 
where the integral is taken over the space occupied by the fluid. We can now utilise condition (|55p to obtain. 



/i = I dS [ hahb'n°'i!'n'V c. ( P + ^ 



I " nb I c 1 m A k 1 c 

+ —UbB rjcb^ - —ecmkn A 



In this expression the first term is symmetric with respect to 77 and ^ while the third term involves only exterior flelds. 
Turning to the second surface term S2 we have, 

l2{v,0 = J dVS2 ^--^ j dSfiaB^r)" {b,+^\bBc) 
Adding both integrals, 

I1 + I2 = I dS hanbr]'^^''n''V c(p+ ^\ ~ -^^crnkn^A'^bl + 



(62) 
(63) 

(64) 
(65) 
(66) 
(67) 
(68) 

(69) 
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ds 



where we have used eqn. ([61 



r}2 

" " a " cr-7 I I 

naUi-r) t, n Vci P + 7— 
8tv 



1 Ilk tc 

47r 



in 



(70) 



Noting that '^VtB^ = Vt-Bc = ^cB^ we can symmetrise the third term, 



vT^bB^c = ^v'Ci^tB^ + ^cBl) 

FinaUy, we can manipulate the 'exterior' second term, 



7 o " cmk A 7 X 

db rim^ AkOc 



dVVr^ 



{e^'^^Akbl) = I dV {hi - Ake 



kmcT-7 tx\ 
VmOc ) 



(71) 



(72) 



Note that the integration in the former two expressions is taken over the exterior vacu um. Since we generally assume 
e"'"^V6-Bc = in the exterior the second term in (I72p vanishes. It it worth noting that Kovetz ( 19661 ) makes the same 
assumption. This is, of course, a simplification since it means that we are not allowing electromagnetic waves. While this 
is consistent with our approximate treatment of the surface, it would obviously be inter esting to relax 1: hese constraints in 
a future analysis. The natural strategy for such a study would be to follow the work of Friedman 1 19781 ') who analyses the 
non-magnetic problem in full general relativity. The magnetic problem ought to be quite similar, with electromagnetic waves 
incorporated essentially in the same way as gravitational waves. 
Collecting the above results, we arrive at 



Im{v,0 = h+h 



dS 



a j-b " ct—7 / 1 ^ 

nanbV ? Vc ( p + — 



Stt 



{haB^)7l't{VbB: + VcBl) 



1 

47r 



dVhi 



and clearly, Im{v,0 = ^m(C,^)- 

Hence we have proven that the remaining magnetic surface term in 



is Hermitian 



(73) 



(74) 



This concludes the proof that the entire C operator in (|22|) is Hermitian in the magnetic problem. It then follows from eqn. (|27|l 
that the symplectic product W{rj,^) is a conserved quantity. 



3 CANONICAL ENERGY AND ANGULAR MOMENTUM FOR MAGNETIC STARS 



As discussed bv lFriedman fc Schuta (|l978al ). the symplectic product W can be used to define the system's canonical energy, 

(75) 

Written explicitly. 



Ecio dv {p\dte + (e")*CHKa]) +/M(r,c) 



(76) 



£2 
47r 



2pi?e[(5")*V6ClKVc«a + V.$} + -L|B°Vaei' + ^B''Re{Va{ey{i!"^cB„ - B^V ^t)] 



B^ 

Stt 

6 



47rG 



dS 



4% 



B^ 

8-K 



2-K 

- ^(rB.)(r)*e'(V5B: + Vcsn 



+ ^ / dVbt 



A conserved canonical angular momentum Jc can be defined in a similar fashion. We have 



(77) 



(78) 



which is identical (as a functional) to the angular momentum expression of the non-magnetic problem ijFriedman fc Schutz 



1978a|). This is not surprising, since Jc does not depend on the C operator. 



It is useful to compare our result for Ec to the corresponding expression for the non-magnetic problem. The latter is 
simply obtained by setting _B' = in the above equations. We see that, in the magnetic case the canonical energy receives 
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contributions from the fluid interior, the vacuum exterior and the fluid-vacuum interface. This makes sense as this energy 
accounts for both fluid and magnetic-field degrees of freedom and the latter extend beyond the fluid body. 

The canonical energy functional plays a central role in the formulation of a stability criterion for a fluid configuration 
coupled to a radiative field. In its original form the criterion states that the configu ration is stable if _Ec(C) > for a/Zinitial data 
( Bernstein et al.lll958l : iLvnden-BeU fc Ostrike^ll967l : IPriedman fc Schut j|l978bl ). In order for the system to be dynamically 
unstable we must have Ec{^) = 0, and if -Ec(C) < for any nonaxisymmetric data, then the configuration is secularly unstable. 
Intuitively, this criterion makes sense since the emitted radiation causes dtEc < invariantly. A perturbation which initially 
has negative canonical energy has the capability of growing in time while further decreasing its canonical energy. From our 
expression for the canonical energy (|77|l we learn that some of the terms associated with the magnetic degrees of freedom 
are positive definite, and thus act to stabilise the system. At the same time, there are magnetic contributions to which 
are indefinite and which may (depending on the field configuration and the displacement vector) stimulate the onset of an 
instability. However, it would be premature to draw any conclusions about stability at this point. The above criterion is 
incomplete. Before applying it we must ensure that the displacements we use correspond to a true physical perturbation. 



4 TRIVIAL DISPLACEMENTS 

An important (and complicating) property of the Lagrangian perturbation t heory for a fluid system is th e existence of a 'gauge' 



freedom associated with the displacement vector As first discussed bv lFriedman &i Schuta (jl978al ). the theory admits a 



class of trivial displacements rj" which correspond to vanishing Eulerian perturbations for all dynamical variables. Although 



the tri vial displacements are unphysical they can have significant impact on the energy stability criterion IjFriedman fc Schutz 



1978a ). This is a consequence of the fact that trivial displacementqf] are associated with a nonvanishing Ec- If we take any 



two displacements and ^2 that differ by a trivial, i.e ^5 = Ci + they obviously describe identical physical perturbations. 
Yet, Ec{£,i + J?) 7^ Ec{£,2) and it is easy to see how this could invalidate any stability conclusions drawn from the canonical 
energy. 

In order to make sense, the stability criterion must be reformulated in terms of 'canonical' displacements ^ which should 
have the property Ec{^ + v) ~ Ec{0- The family of canonicals is formally singled out by the requirement of orthogonality 
with respect to the product W , 

W{i, T)) = 0, for all trivials 77 (79) 

It follows that if the canonicals ^1 and ^2 describe the same physical perturbation then ?7 = 1^2 — ^1 is a trivial and since both 
^1 and ^2 are orthogonal to all trivials, any trivial is itself orthogonal to all trivials, W{ri, rj) = 0. Then clearly, 

£c(6) = i w^(a6, 6) = I widtiu^i) + \ w{dtn, v) + w{dtii,v) = ^ w(dtii,ii) = i5c(6) (so) 

Now that we understand the nature of the problem, it is clear that we need to investigate the class of trivial displacements 
permitted in the magnetic problem. 

The general definition (|79|l leads to an explicit condition satisfied by every canonical displacement. In order to derive 
this condition we first need a solution for the trivials. In the fluid problem, the trivial diplacements 77" are defined by the 
requirement Sp = Su"" — which means, 

Vaipvl = (81) 
[dt + r^Jr?" = (82) 

When a magnetic held is present we also require &° = which provides one additional equation 

t^'^Vi, {e^aeV'^B') = (83) 

The remaining perturbations 5p and 6E°' vanish as a consequence of the equation of state p = p{p) and eqn. ([5]), respectively. 
Note that according to (|82p trivial displacements are always Lie-transported by the background flow. In fact, the relation (I82II 
can be used to show that 

W{i v)= I dVpr)" A,Ua = (84) 



Before moving on to the m agnetic problem, it is u seful to recall what happens for a non-magnetic fluid. Ignoring (|83l 
the solution of (dU and ([82| is (jFriedman fc SchutzllmSai ) 



^ Note that for any trivial displacement -Bc(»?) = 0, provided the background configuration is static (see Appendix A of lFriedman &: Schutd 
In this case there is no issue of gauge dependence for the canonical energy. 
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jf = le^'^VbCc (85) 

p 

where the vector field C,c is arbitrary, as long as it is Lie-transported along the background flow, 

[dt + jCuKa = (86) 
This last result follows directly from 5u'^ = and the identity 

CuiVaFt) = VaiC^Ft) - F.VaVbu" (87) 

Using (f85)) in ((84)) we obtain, 

dVe'''"'Vb{A^Ua)(c = (88) 
Since C^c is arbitrary, the above condition is satisfied if 

e'""^Vb(A|Me) = Aj (^^e'""'VbUc^ = (89) 

The vector appearing in the last expression is the fiuid's specific vorticity uj"'/p = e°'"^S7bUc/ p. Originally derived by 



Friedman fc Schuta (|l978al '). expression (|89|l is the general condition which should be satisfied by all canonical displacements. 
It states that a canonical displacement must preserve the vorticity of each fluid element. 

A different aspect of the same result can be revealed by taking the curl of the non- magnetic Euler equation ([T]). This 
leads to the familiar vorticity equation, 

— UJ Vbu"' — Lj'^VbU (90) 

This is easily written as, 

[a + (^y ^ = (91) 

which shows that the specific vorticity vector is Lie-transported. Perturbing (|9ip . 

[at + £.u]A (^y^ =0 (92) 

Combining this general result with (|89p we see that the canonical displacements single out the the trivial solution A(lj"/ p) = 
of eqn. (|92|l . 

Finally, it is interesting to note that eqn. (|92|l can also be derived from dtW{S,,rj) — (where now ^ is not necessarily 
canonical). Starting from (|84p we can show that, 

W{^,7l) = J dVCcpA(^j^ (93) 
Taking the time derivative of this expression and making use of (|86l) we find 

dtW{^, n)= j dVCcP [dt + A,] A (^^) = (94) 

Since is arbitrary we end up with eqn. (|92|) . 

Having arrived at the condition for a canonical displacement, we still face the issue of how to determine such displacements 
in practice. For later convenience, we note that 

A^Ua = Vaf (95) 

where / an arbitrary scalar function, is always canonical 0. Indeed, 

VK[e,r?]=0 ^ [dVVa{pvlf = (96) 



is identically satisfied due to (|8ip . 

Moving on to the magnetic problem, let us first provide an indication that it may be somewhat challenging. It is easy 
to see that the trivial form (|85|) with a general ("^ cannot possibly be valid in the magnetic case. Indeed, if it were to remain 
true we could reproduce the derivation of (I92p . However, in the MHD framework eqn. (I9ip is replaced by a more complicated 
relation, 



^ This can be immediately verified from l|89|l for a non-magnetic fluid, and the same is true for the case of a non-isentropic fluid 
llFriedman fc Schut jri978ah . 
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+ (^)^^M'' (97) 
with 

= e'""^ B^VtEk Vc(l/p) + e^'^^Vb (^i 5*= V^Bc j (98) 

Since the right-hand side of this equation is non- vanishing, we see that the specific vorticity is not conserved in MHD. In fact, 
as far as we are aware, no corresponding conservation law is known. It is easy to see that this may complicate the analysis of 
the canonical displacements. 

The equations that need to be solved in order to identify the class of trivial displacements are (|81|l and (I83|l . Since the 
general solution to the first of these equations is already known from the fluid problem we 'only' need to solve the magnetic 
equation and then link the two solutions in the appropriate way. First note that the general solution to (|83|l can be written 



e,aeV''B^ = Vci/- (99) 
where Tp is an arbitrary scalar functioij^. This equation shows that 

r^'^Vai^ = B'Va^ = (100) 

In other words, both the trivial displacement and the magnetic field lie in surfaces of constant -0. Taking the cross product of 
((99)) with B" we obtain 

5^77" - (r?*B5)B'" = e'''"'BtVc^ (101) 

Given this, it is easy to show that if we decompose the trivial into r;" ~ + Vli where rj^^ and r;" are respectively parallel 
and orthogonal to B", then 

r?i = ij,-^-B,V,i^ (102) 
We cannot, however, say anything about the parallel component at this point. Representing it by a scalar function a, we have 
77" = aB" + -ije'""=Bi,Vc^ (103) 
This expression must also satisfy the fiuid trivial condition (|8ip . This leads to a relation between a and tp, 
B'^Vaipa) + e^'^Va (^) (VcV) = (104) 



V B2 . 

Unfortunately, this relation shows that one should not expect to be able to completely determine a from a given if) (or vice 
versa). To see this note that the first term in (|104|l is a directional derivative along B\ while the second term is a scalar 
function (let us call it g). Suppose that we introduce a coordinate system [61,62,63] such that e\ is aligned with B\ Then 
(11041) can be written 



B^^+g{xi,X2,X3)^0 (105) 
which integrates to 
„ _ 1 /■ 5 



dxi + h{x2,X3) (106) 

P J B 

where the function h is freely specifiable. From this simple demonstration we learn that we need to treat the functions a and 
ip as (at least partially) independent. 

A simple example of a case where a and ip are independent is a uniform density fluid and a uniform force-free magnetic 
field. Since e'''"'VbBc — for a force-free field and taking p and B^ as constants, the constraint (|104[) collapses to 



B^Vcipa) = (107) 

from which it is clear that ip remains freely specifiable. 

Even though we have failed to obtain an explicit expression for the trivial displacements in terms of a single arbitrary 
function, we can proceed to write down the condition that determines a canonical displacement. Following the steps that led 
to (1891) we arrive at 



^ In terms of the magnetic vector potential A we have 
b = V X 5A = V X (5 X B) 

which gives SA = g X B + V/ where / represents the gauge degree of freedom. For a trivial displacement we have r; X B = Vi/i hence 
5A = V(/ -I- Hence, the trivials in the magnetic case generate a pure gauge transformation for the magnetic vector potential. 
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paB" - Vc 



V B 



CO 



P 



A^ua - V j A^- ( — ) = (108) 



Although this is not quite satisfactory, when combined with a solution to pU4|) it does allow us to confirm whether a proposed 
perturbation is canonical or not. In fact, we already know that a class of canonical perturbations are given by (|95p . It is easy 
to show that this class remains acceptable also in the magnetic case. 

As an example of applying (|108p let us work out another set of canonical displacements. We see from (|108|l that the 
problem simplifies if we take 

B^A^Ua = (109) 
To satisfy (|108|l we additionally need (since is general) 

(^^0) 

One can show that this is satisfied by 

A^Ua^-iB^VaXc-B^VcXa) (HI) 
^ P 

for a general Xa- 



5 STABILITY CONSIDERATIONS 



Having at hand the canonical energy functional -Bc(5) and a properly formulated cri terion in terms of canonical diplacements, 
we are ready to investigate the stability of various model systems. As discussed by Friedman fc Schut j ( 1978ah this analysis 
can be based on either normal modes or canonical initial data. In the case of normal modes, i.e. solutions taking the form 
^ ~ gSCTt+im^ (with m the usual integer eigenvalue associated with the separation of the azimuthal (p coordinate and a the 
frequency) one can show that they generally correspond to canonica l displacements. Although there are exceptions to this 
rulqj an instability proof based on a mode-solution is generally valid ( Friedman &: Schutall978al ). 

Another important result that remains valid in the magnetic problem is the relation between the canonical energy and 
angular momentum 

^=-P = -- (112) 

This shows that the onset of instability can be deduced from a change in sign of the mode pattern speed Up. If we imagine 
the star spinning up, an initially counter-rotating mode with Jc < becomes unstable (in the sense of having negative Ec) 
beyond a critical rotation frequency. The proof of the above relation proceeds exactly as in the non-magnetic problem since 
Ec and Jc have the same functional form when expressed in terms of the A and B operators of the master eqn. (|22|l . What 
is different between the magnetic and non-magnetic problems is the intrinsic character of the modes themselves. A stability 
analysis for various classes of oscillation modes requires the determination of the associated displacement vectors. Solving 
for the oscillation modes of a rotating magnetised star is, however, a challenging problem. Once this problem is solved, the 
canonical energy formula can readily be used to test stability. 

In absence of actual mode-solutions for a rotating magnetic star, we can make progress by discussing various classes of 
initial data. This will provide some insight into the role of the magnetic field. Let us start by extending the demonstration 
that a rotating star is generically unstable due to the emission of radiatio n . In th eir discussion of the secular instability 
of (isentropic) fiuids coupled to gravitational radiation, Friedman fc Schutz ( 1978bh chose the following canonical data (in 
cylidrical coordinates to, z, <p), 



t = {pS; + !.")e""*' 
with / = Ae""*'. Given this choice we find. 



A^Ua 



V/-a i: 



U Vbt, 



(113) 
(114) 

(115) 
(116) 
(117) 



^ In the case of a fluid configuration with purely azimuthal flow u''' = a normal mode can belong to the family of trivials if it is 

corotating with the background fiow, i.e a = —mQ. 
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\dte 



1 



AA* + {u^Y{'iuj^ny* + XX* + ZZ*) - 2M^(/iA* + ^*A) 



(118) 
(119) 



where Z = and X = . The resulting expression for the perfect fluid canonical energy, minimised with respect to A, is 
Ec < —m I dVpyCg ~ uj Q j/j,^ + mfci + fc2 



(120) 



where — ^p/p is the sound speed and the functions fci and ^2 are independent of m. Since Cs — > near the surface of the 
star, the energy functional becomes generically negative in the 'multiarm' limit, m S> 1. Therefore, we can find physically 
acceptable perturbations that make all rotating fluid stars unstable to gravitational radiation. 

We can use the same initial data, eqns. (|113p & (|114p . in the full MHD canonical energy ([77}. In addition to (|115p - H119[l 
we then need the magnetic terms. 



m\B'^f{zu'^fifi* + XX* + ZZ*) 

m^B*' {ccjVm'-B^ + Re [m(S"X + B"" Z)] } 



(121) 
(122) 



This leads to 

1 
2 



Ec < \rn I dVp \ pp* { c, - m^n"^ + ^ ) + 



{Bn 



47rp I Anp 



■(XX* +ZZ*) 



^ [B^{p*X + pX*)+ B''{p*Z + pZ*) ] + mfci + fca 
47rp ' 



(123) 



where = {B^Y + {B^Y'. Since the first two magnetic terms are positive definite they always play a stabilising role. 
Meanwhile, the third magnetic term can be either positive or negative. Hence, it could potentially destabilise the system. 
It should be noted that this term vanishes when the magnetic field is either purely poloidal [B'^ = 0) or purely toroidal 
(B^ = B' = 0). It should, however, be present for all realistic neutron star field configurations. In general the (negative) 
rotation term is the main destabilising factor. The main stabilising term is 7p/p which is a surrogate for the star's gravitational 
potential energy. This energy is far greater than the energy stored in the magnetic field and, consequently, provided we are 
not too close to the surface, the sound speed term will dominate over any magnetic term (assuming components for the 
displacement vector p,X, Z of roughly the same order of magnitude). Hence, one would expect the instability region to shift 
only slightly in the magnetic case. If a non-magnetic star is unstable at a given rotation rate for m S> 1 then it is very likely 
that the same is true for the corresponding magnetic star. 

To conclude t he discussion, let us consider a specific case in spired by the instability of the inertial r-modes in a rotating 
perfect fluid star (Andersson 19981 : Friedman fc Morsink 1998f ). Inertial modes are incompressible to leading order in the 
slow-rotation expansion. This means that it is interesting to consider displacements for which Vi^' = 0. In this case, the 
canonical energy simplifies considerably. Let us, for the sake of the argument, assume an incompressible mode-solution of 
form (|113|) with time-dependence represented by a frequency o (although it should be pointed out that it is not at all clear 
that such a solution actually exists). Then, further restricting the analysis to a purely azimuthal magnetic field, we find in 
the m ^ 1 limit. 



Ec < 



dVp 



o 

pp + V Vt) 



[B 



Anp 



+ mki + k2 



(124) 



This result is interesting because it suggests that in order for the mode to have a chance of becoming unstable the star's 
rotational frequency must first dominate over the Alfven frequency of the background magnetic field. This is a strong hint 
that the instability of (say) the r-modes may be s ignificantly affected by the presence of a magnetic field. This adds further 



_Dy tn 

complexity to the interesting issues discussed bv iRezzoUa et al.l l|2001alla) . The problem of radiation driven instabilities in 



magnetic stars undoubtedly warrants further investigation. 



6 CONCLUDING REMARKS 

In this paper we have presented a framework for analysing the stability properties of a stationary and axisymmetric, infinitely 
conducting perfect fluid conflguration threaded by a magnetic held and surrounded by vacuum. We have developed a La- 
grangian perturbation approach which enables us to formulate a rigorous stability criterion based on the canonical energy. 
For any choice of initial data, represented by the displacement vector and its time derivative {£,, dt^}, the sign of the canonical 
energy determines whether the configuration is stable or not at the linear perturbation level. 

Our analysis provides the first complete description of the stability problem for a rotating magnetic star. We extended 
previous results for rotating stars by allowing for the presence of an exterior vacuum field. We also discussed (for the first 
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time in this context) the comphcations associated with the so-called 'trivial' fluid displacements, which do not represent true 
physical perturbations. In order for the stability criterion to make sense one has to isolate these trivials and consider only 
the physical 'canonical' displacements. Wc discussed this problem and formulated a condition which must bo satisfied by all 
canonical displacements. Having obtained a well-defined stability criterion we provided two examples that indicate that the 
magnetic field tends to stabilise radiation driven instabilities. 

To make further progress we need to calculate explicit mode-solutions for a rotating magnetic system. Once such modes 
have been obtained the canonical energy formulated in this paper can be used to assess their stability. Of course, the calculation 
of oscillation modes of a magnetic star is a far from trivial problem that sets a severe challenge for future efforts in this area. 
On the other hand, strong motivation for such work is provided by increasingly detailed observations of for example magnetar 
flares and associated quasiperiodic oscillations. 



APPENDIX A: PROOF OF HERMITICITY OF THE C OPERATOR 

In this Appendix we describe the various steps involved in proving the Hermiticity of the C operator. We start from (see the 
main text) 



47r 



or. 



■^b'Be 
An 



-VaTl'" [dp+^b'Bc^ - ^n^b^VcBa + B'V ,ba) " P»7TVc(m'' VfcMa) 



Our task is to manipulate ?;"C[^a] to make it symmetric with respect to rf and 
Let us first study separately the magnetic 'tension' term, 

- ^rt^ib^VcBa + B'Vcba) = -^V^i^'^^bB^VcBa - B'-VbC^cBa + VbtB^VcBa) - -^V civ" B^ba) 

+ l-B'VcV"i-^b£.''Ba - ^\bBa + B^Vb^a) 



(Al) 



(A2) 



-^Vc{v''B%a) + ^B'^B'^Vcff'^Vbia - -^VbtB^iBaVrV " V'^^cBa 



4tv 



1 

An 



47r 



{-Tft^bB^VoBa + ri'^B^be^cBa + B"Vc»?"$'V6B„) 



And since, 

jfB^VbC^cBa = Vb{v''B''eVcBa) - VbrfB^^VcBa - rfCB^bVcBa 
we have. 



47r 



rfib'VcBa + B^Vcba) 



~Vc{v^B%,+n^B^^''VbBa) + ^B^B'V.n^Vb^a 



VbS^B^iBaVc-n" - V'^cBa) + ^ri''C^b{B''V.Ba) 



1 

47r Ait 
Using once the background Euler equation in the last term, wc arrive at 



4lT 



rj^ib^WcBa + B"Vc6a) 



-^Vc{v''B%a+v"B^e^bBa) + ^B'B^VcV'Vb^a 
iTT 47r 



1 



-— Vi,fS=(B„Vcr?" - 7?"VeBa) + ?7"f V„V6 [P+-^ 



B' 



+»?"rV6[p(w=VcWa+Va*)] 

Moreover, we have for the total pressure term, 

- Va??" {Sp + ^b^Bc^ = VaV'ilP^bt + t^bP) - ^VaV^Bci-Vb^'B' - ^"VbB" + B'VbC 



"^aV'^^bti-yP + B^An) + WaV't'^bip + B^Sn) - -^S/ an'' B^V bC Be 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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Using once more the background Euler equation in the second term we obtain, 
Collecting the above results and returning to the full expression, 



(A8) 



47r 



47r 



+v"^''^aVbip + b^/8tv) + Vav'^bti'yp + b^/^-k) + pn'-VaS^ 

-pV^e'^ciuVbUa) + pn'^iu^X/bfia + Jy^C" Vi, [p(^i'= Velta + V„$)] - pVaT/^^ ' ( Vfc* + u'V ,Ub) (A9) 

The remaining non-symmetric terms (apart from the first surface term) are the last five terms. Denoting them collectively as 
R we have, 

R = -r;"V6(pC'')(V„$ + u^V ^Ua) ~ pv^C'^oiuVbUa) + pv'"iu\7i)\a 

Wt^bHuVcUa + V,*)] - pVaV'^^'^b^ + uV^Ut) (AlO) 

Differentiating by parts the fourth term we find, 

R = -piu^V.Ua + V„$)(r?''V6e' + C^bV") + PV't^aVb^ + pV^uVbiuV^C) (All) 

Furthermore, 

prfu'Vbiu'V cC') = p'n°'u'ybu''Vc^a+pr)°'uu''Vbycia 



= VbiP'n^Uu'V c(,a) — pV bV" U u'V da ~ daiuV bP + pV bu) 

The last term vanishes due to the background continuity equation so eventually we are left with, 

R = Vi,(p77%'«'V4a) - p{u\I.Ua + Va $) (?7° Vt^" + TViT?") 



(A12) 



(A13) 



Written in this form, 7? contains only su rface and symmetric terms. The final step is to symmetrise the gravitational potential 
term prj'^VaS^ in eqn. (|A9|I . From Ref. [Friedman fc Schutj (11978a|) we have. 



V47rG 



^5^$ + pr," 



(A14) 



Collecting all the results we finally obtain. 

1 

4tt 



ri''{Sp + B%c/4tt) 



V^iB^, + B'^^'VbBc) + pri^^u'V^b + (4^^"^"* + 



1 

'4^ 



V^Sr^^Vbh^ + ^B^B^y^Tj^Vbia + ??"C''V„V6(p + B'/Stt) + V b^ {iV + BV47r) 
47r 



47r 

-p(^i'=VcU, + Va$)(77'"VbC'' + C^brt) + pr?"C'V„ Vi,$ - pu"ii'Vi,r?" V4a 
This is eqn. (|28|) of the main text. 



(A15) 
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